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In this work we compute models for relativistic white dwarfs in the presence of strong magnetic 
fields. These models possibly contribute to super-luminous SNIa. With an assumed axi-symmetric 
and poloidal magnetic field, we study the possibility of existence of super-Chandrasekhar magnetized 
white dwarfs by solving numerically the Einstein-Maxwell equations, by means of a pseudo-spectral 
method. We obtain a self-consistent rotating and non-rotating magnetized white dwarf models. 
According to our results, a maximum mass for a static magnetized white dwarf is 2.13 Mq in the 
Newtonian case and 2.09 Mq while taking into account general relativistic effects. Furthermore, we 
present results for rotating magnetized white dwarfs. The maximum magnetic field strength reached 
at the center of white dwarfs is of the order of 10^® G in the static case, whereas for magnetized 
white dwarfs, rotating with the Keplerian angular velocity, is of the order of lO'^"^ G. 


I. INTRODUCTION 

White dwarfs (WDs) are stellar remnants of stars with 
masses of up to several solar masses. With a mass compa¬ 
rable to that of the Sun ~ 2 x 10^^ g, which is distributed 
in a volume comparable to that of the Earth, the central 
mass density in these objects can reach values of about 
10^^ g/cvo?. Together with neutron stars and black holes, 
they are the endpoints of stellar evolutions and play a key 
role in astrophysics 0,i- 

The existence of white dwarfs was one of the major 
puzzles in astrophysics until Fowler Q, based on the 
quantum-statistical theory developed by Fermi and Dirac 
[j, Q , shows that white dwarfs are supported by the pres¬ 
sure of a degenerate electron gas. In addition, Chan¬ 
drasekhar in Ref. @ included effects of special relativity 
in the degenerate electron gas theory and, as a result, 
found that there is a limit in the stellar mass, above which 
degenerate white dwarfs are unstable. This critical mass 
is the so-called Chandrasekhar limit and is about 1.4 Mq. 

White dwarfs are mostly composed of electron- 
degenerate matter. The mass of the star is essentially 
due to the nuclei, whereas the main contribution to the 
pressure comes from the electrons. For typical WDs, 
thermonuclear reactions terminate at lighter nuclei, as 
Carbon, helium, or oxygen d, Q . 

Some white dwarfs are also associated with strong 
magnetic fields. From observations, the surface magnetic 
field of these stars can reach values from 10® G to 10® G 
[sl-fT^. However, the internal magnetic field in magnetic 
stars is very poorly constrained by the observations and 
can be much stronger than in the surface. For example, 
white dwarfs can have internal magnetic fields as large as 
2 ^q 12-16 q according to Refs. [U, [ij, [3 ■ Moreover, self- 
consistently highly magnetized neutron stars calculations 
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have shown that neutron stars can possess central mag- 
netic fields as large as 10^® G |15l - [la| . Therefore, under¬ 
standing and estimating magnetic fields inside compact 
objects occupy a key position in astrophysics. 

Motivated by observations of a thermonuclear super¬ 
nova that appears to be more luminous than expected 
(e.g. SN 2003fg, SN 2006gz, SN 2007if, SN 2009dc), 
it has been argued [l9l - l^ that the progenitor of such 
super-novae should be a white dwarf with mass above 
the well-known Chandrasekhar limit, in other words, a 
super-Chandrasekhar white dwarf. 

Progenitors with masses M > 2.0 Mq were considered 
in the literature as a result of mergers of two massive 
white dwarfs, or due to fast rotation 3- In addition, 
super-Chandrasekhar white dwarfs were investi gat ed in 
a strong magnetic field regime as in Refs. [25l-l3- In 
the Newtonian framework, models for white dwarfs en¬ 
dowed with a magnetic field and/or rotating were inves¬ 
tigated in a series of papers [2814^ . Recently, a study 
of differentially rotating and magnetized white dwarfs 
were performed within the ideal magnetohydrodynamic 
(GRMHD) regime [3 nnd they have found that differ¬ 
ential rotation can increase the mass of magnetized white 
dwarfs up to 3.1 Mq. 

In Ref. [3 j it was obtained ultramagnetized white 
dwarfs with a maximum mass of 2.58 Mq. The au¬ 
thors also found a magnetic field strength up to 10^® G 
at the center of the star. Nonetheless, such approach vi¬ 
olates not only macro physics aspects, as for example, 
the breaking of spherical symmetry due to the magnetic 
field, but also micro physics considerations, which are rel¬ 
evant for a self-consistent calculations of the structure of 
these objects Besides, a self-consistently newto- 

nian structure calculation of strongly magnetized white 
dwarfs shown that these stars exceed the traditional 
Ghandrasekhar mass limit significantly (M ^1.9 M«) for 
a maximum field strength of the order of 10^“^ G [TJ. 

In this work, we model static and rotating magnetized 
white dwarfs in a self-consistent way by solving Einstein- 
Maxwell equations in the same way as done originally for 
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neutron stars in Refs. We follow also our recent 

work on highly magnetized hybrid stars [l^ j where both 
general relativistic effects and the anisotropy of the en¬ 
ergy momentum tensor caused by the magnetic field were 
taken into consideration to calculate the star structure. 
The presence of such a strong magnetic field can locally 
affect the microphysics of the equation of state (EoS), 
as for example, due to Landau quantization. As shown 
in Ref. [l^, the Landau quantization does not affect the 
global properties of white dwarfs. However, the authors 
solve the structure equation in a Newtonian form and 
without the magnetization term for the matter. More¬ 
over, as we will see, effects of general relativity can play 
an important role by determining the maximum mass of 
such highly magnetized white dwarfs. 

Globally, the magnetic field can affect the structure 
of WDs since it contributes to the Lorenz force, which 
acts against gravity. In addition, it contributes also to 
the structure of the spacetime, since the magnetic field 
is now a source for the gravitational field through the 
Maxwell energy-momentum tensor. As we are interested 
in global effects that magnetic fields and the rotation 
can induce in WDs, we simplify the discussion assuming 
white dwarfs composed predominately by {AjZ = 2) 
in a electron background . 


II. BASIC EQUATIONS AND FORMALISM 


In this work, we consider rotating and non-rotating 
magnetized white dwarfs. The formalism used here was 
first applied to stationary neutron stars [H, [H, [s^ and 
more recently in Ref. [^. Details of the gravitational 
equations, numerical procedure and other properties of 
the equations can be found in the references cited above 
and in Ref. [s^- For the sake of completeness and 
for a better understanding of the reader, we show here 
some of the electromagnetic equations that, together with 
the gravitational equations, are solved numerically. The 
energy-momentum tensor of the system reads: 

Tap = {e+p)UaUp+pgaP + — (Fa^F!^ “ ^F^^F^^gap] , 

po \ ^ 4, J 

( 1 ) 

where Fap, denotes the antisymmetric Faraday tensor, e 
the energy density and p the pressure as measured by an 
observer (Oi) co-moving with the fluid and whose four- 
velocity is Ua- The gap is the metric tensor. The first 
term in Eq. ([IJ is the isotropic matter contribution to 
the energy momentum-tensor, while the second term is 
the anisotropic electromagnetic field contribution. 

The electromagnetic contribution (EM) to the energy- 
momentum tensor is obtained within the so-called 3-1-1 
decomposition [H, The energy density becomes: 

+ ( 2 ) 

2fio 


and the momentum-density flux can be written as: 

j{EM) ^ ^3^ 

^ Po 

The stress 3-tensor components are given by: 

S(EM)r = J_(^E^Ee - E’'Er + B^Be - B^Br), (4) 
2^0 

g{EM)e ^ 1 _ B^Be), ( 5 ) 

2po 

giEM)p ^ ^^E^E, + B^B,), ( 6 ) 

being the electric field components, as measured by the 
Eulerian observer Oq , written as [38j| : 
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and the magnetic field given by: 
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( 7 ) 
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with N^{r,B) being the shift vector, N{r,B) the lapse 
function and T a metric potential (for more details see 
Refs, As in Ref. , the equation of motion 

= 0) reads: 

E[ (r, B) + 1 / (r, B) — InT (r, B) + M (r, B) = const, (9) 


with F[ (r, B) being the logarithm of the dimensionless rel¬ 
ativistic enthalpy per baryon: 


H := In 


e+P \ 
nihnbc^ ) ’ 


( 10 ) 


with mb the mean baryon mass 1.66 x 10 ^^kg, nb the 
baryon number density. The second term in Eq. is 
defined as i/ = v{r,B) := ln{N) and the Lorenz factor 
written as T = (1 — [/^)“ 2 . The physical fluid velocity 
U in the (j) direction is defined as: 


u = (11, 

and the magnetic potential M(r,B) associated to the 
Lorentz force is written as: 

M (r, B)=M {A^ (r, B)) ■.= - f f {x) dx, (12) 

J A,l,{r,9) 

with a current function f{x) as defined in Ref. (see 
eq. (5.29)). In this work, we use a current function 
f{x) = fo = const, which is proportional to the intensity 
of the magnetic field, i.e., the higher the current function 
/o, the higher is the magnetic field in the star. As shown 
in Ref. [i3, other choices are possible for f{x), however, 
they do not alter the conclusions. 
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III. MASS-RADIUS DIAGRAM FOR STATIC 
HIGHLY MAGNETIZED WHITE DWARES 

In this section we present the mass-radius (MR) di¬ 
agram for static magnetized white dwarfs. The rela¬ 
tion between the mass and the radius for non-magnetized 
white dwarfs was first determined by Chandrasekhar Q . 
Recently, studies on modified mass-radius relation of 
magnetic white dwarfs were proposed, for example, in 
Refs. [l3, [H, [13. As we also found in this work, those 
authors shown that the mass of white dwarfs increase in 
presence of magnetic fields. 

In Fig. m we show the isocontours in the {x, z) plane 
of the poloidal magnetic field lines for a static star with 
central enthalpy of He — 0.0063 c^. As we will see in 
Fig. m this value of the enthalpy gives us the maxi¬ 
mum gravitational mass of relativistic, static and mag¬ 
netized white dwarfs achieved with the code, namely, 
2.09 Mq, which corresponds to a central mass density 
of 2.79x 10^° g/cm^. It is known that at sufficiently high 
densities reactions as inverse /3-decay or pycnonuclear fu¬ 
sion can take place in the interior of white dwarfs [13,113 ■ 
As estimated in Ref. [13, for the maximum mass config¬ 
uration obtained in this work, i.e., for a central magnetic 
field of 1.03 x 10^^ G, the onset of electron capture by 
carbon-12 nuclei was found to be about 4.2 x 10^° g/cve? 
with electron-ion interactions, and 3.9 x 10^° g/cm^ with¬ 
out electron-ion interactions. Therefore, in our calcula¬ 
tion, the maximum mass density reached by the most 
massive and non-rotating magnetized white dwarf lies 
below the threshold density for the onset of electron cap¬ 
tures by carbon-12 nuclei as calculated in Ref. [13 • 

In Fig. [2 we show the mass density distribution for 
the same star as shown in Fig. [TJ As expected, the mass 
density is not spherically distributed and the maximum 
mass density is not at the center of the star. In this case, 
the central magnetic field reaches a value of 1.03x 10^® G, 
whereas the surface magnetic field was found to be 2.02 
XG. As a result, one sees that the Lorentz force ex¬ 
erted by the magnetic field breaks the spherical symme¬ 
try of the star considerably and behaves as a centrifugal 
force, which pushes the matter off-center. In order to 
understand better this aspect, we can use the equation 
of motion m for the static case F = 0: 

H{r,9) + + M{r,9) = C, (13) 

and then plot these quantities in the equatorial plane as 
shown in Fig. [3) This same analysis was presented for 
magnetized neutron stars in Ref. Ea. The constant C 
can be calculated at every point in the star. We have 
chosen the center, since the central values of the mag¬ 
netic potential M{r,9) is zero and the central enthalpy 
He is our input to construct the models. The Lorentz 
force is the derivative of the magnetic potential M{r,9) 
in the equation Eq. (1131) and reaches its maximum value 
off-center {veq 350 km, see Fig. [S]). As already dis¬ 
cussed in Ref. Ea , the direction of the magnetic forces in 
the equatorial plane depends on the current distribution 


Magnetic field 



X [km] 


FIG. 1. Isocontours of the magnetic field strength in the {x, z) 
plane, with a gravitational mass of 2.09 M© and a magnetic 
dipole moment of 1.30 x 10^“^ Am^. The ratio between the 
magnetic pressure and the matter pressure at the center of 
the star is about 1, and the magnetic field at center reaches 
1.03x10^® G. 


baryon number density 



-500 0 500 


X [km] 

FIG. 2. Isocontours of the baryon number density 
in the {x, z) plane for the same star as shown in 
Fig. m The central baryon density for this model is 
1.679 X 10“® fm“® (2.79 x 10^° g/cm®). 
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FIG. 3. Behaviour of the different terms of the equation of 
motion as a function of the equatorial coordinate radius for 
the same star as shown in Fig. [T] 


inside the star. In addition, the magnetic field changes 
its direction in the equatorial plane and, therefore, the 
Lorenz force reverses the direction inside the star. In our 
case, this can be seen from the qualitatively change in 
behaviour of the function M{r,9) around req 350 km 
in Fig. [3l 

In order to compare our results with those in the lit¬ 
erature, we compute the mass-radius diagram for mag¬ 
netized white dwarfs as shown in Fig. |4l As pointed out 
in Ref. (4ll| . a fully consistent equilibrium configuration 
for magnetized white dwarfs is still lacking. Those same 
authors have used a relativistic framework with cylin¬ 
drical coordinates and the splitting of the pressure into 
parallel and perpendicular components to show that the 
maximum magnetic field inside these objects can not ex¬ 
ceed 1.5 X 10^^ G. Their solutions also indicate that it is 
not possible to have stable magnetized WDs with super- 
Chandrasekhar masses. The maximum magnetic field 
strength obtained in Ref [4l| is less than the value ob¬ 
tained by a self-consistent solution of Newtonian white 
dwarfs as presented in Ref. [ijfand orders of magni¬ 
tude less th an p redicted in Ref |^ . However, as shown 
in Refs, the work [2^ has been criticized above 

all because their calculation violate both macro/micro 
physics properties essential for the stability of these ob¬ 
jects. 

In Fig. m the maximum white dwarf masses are ob¬ 
tained when the ratio between the magnetic and the mat¬ 
ter pressure at the center of the star is less than or about 
1. For such a strong magnetic field, the magnetic force 
has pushed the matter off-center and a topological change 
to a toroidal configuration can take place [i3- This gives 
a limit for the magnetic field strength that can be com¬ 
puted within this approach, since our current numerical 
tools do not enable us to handle toroidal configuration. 

The maximum masses of white dwarfs increase with 



FIG. 4. Mass-radius diagram for magnetized white dwarfs. 
Different curves represent different values of the current func¬ 
tion /q. We also compare the maximum white dwarf mass 
obtained in the Newtonian case (in black) and in the rela¬ 
tivistic one (in red) for the maximum electric current value 
for which numerical convergence is achieved. This diagram 
is quite similar to the MR diagram calculated in Ref. El- 
However, those authors have evaluated only Newtonian white 
dwarfs. All curves in this figure were calculated for a ratio 
between the magnetic and the matter pressure less than or 
quite close to 1 at the center of the star. 


the magnetic field. In our calculation, we found a rela¬ 
tivistic white dwarf mass of 2.09 Mq for almost the same 
magnetic field strength at the center, B ~ 10^“^ G, as 
in Ref M- The same authors presented configurations 
with magnetic fields up to 10^® G, which we have not 
found in our calculations. In the Newtonian case, we 
found a mass of 2.13 Mq for the most massive magne¬ 
tized white dwarf In both cases, the masses are well 
above the Chandrasekhar limit of 1.4 Mq. 


IV. ROTATING MAGNETIZED WHITE 
DWARFS 

Aside magnetic fields, rotation is a crucial observable 
in stellar astrophysics. In the case of neutron stars, for 
example, the magnetic dipole model states a direct rela¬ 
tion between the rotation period/period derivative with 
the polar magnetic field of the star. For white dwarfs the 
interest in understanding these objects has been increas¬ 
ing with the years. Some observed white dwarfs rotate 
with periods of days or even years. One of the fastest 
observed WD possesses a spin period of 13.2 s [d^, a 
value similar to those observed in Soft Gamma Repeaters 
(SGRs) and Anomalous X-ray pulsars (AXPs), known 
as magnetars [H, |4^ . A relation between white dwarfs 
and magnetars was addressed in Ref. [1^, where the au¬ 
thors speculated that SGRs and AXPs with low magnetic 
field on the surface might be rotating magnetized white 
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FIG. 5. Mass-radius diagram for static and rotating white 
dwarfs. The TOV solution is shown in black, and in red 
(dashed) we show the Keplerian sequence for rotating WDs. 
On the bottom left-hand corner we show the gravitational 
mass as a function of the central density for the same sequence 
of stars. 

dwarfs. 

Rigidly rotating non-magnetized white dwarfs were al¬ 
ready studied long time ago in the Newtoninan frame¬ 
work as in Refs. (46l - l^ . In addition, the structure of 
rapidly rotating white dwarfs was performed in general 
relativity as in Ref. [l^, and more recently in Ref. 
where the authors used the Hartle’s formalism to 
solve the Einstein equations. From the standpoint of ro¬ 
tation, it is clear that all rotating stars have to satisfy 
the mass-shedding, or Keplerian limit, as a condition of 
stability. This limit is reached when the centrifugal force 
due to the rotation does not balance gravity anymore and 
the star starts to lose particles from the equator, defin¬ 
ing an upper limit to the angular velocity of uniformly 
rotating stars. 

In the same way that rotation provides a natural limit 
for the stability of stars, in the case of WDs, there are 
also microphysics aspects, as for example, the inverse 
/3—decay and pycnonuclear fusion reactions [^. l40l|. 

that need to be taken into account for a complete and 
self-consistent description of these stars. In this paper, 
however, we restrict ourselves to the study of the com¬ 
bined effect of rotation and magnetic fields on the global 
structure of WDs and we do not address here the micro¬ 
physics aspects. 

In Fig. we show the Tolman-Oppenheimer-Volkoff 
(TOV) solution for the structure of a spherically symmet¬ 
ric white dwarf and the mass-shedding frequency limit. 
The centrifugal force exerted by the rotation acts against 
gravity, which allows the star to support higher masses 
compared to the static case. In the first place, by com¬ 
parison of Fig. 0] and Fig. [51 one sees that magnetic fields 
are more efficient than rotation in increasing the maxi¬ 
mum mass of stars. The maximum mass obtained for a 


FIG. 6. Keplerian frequency as a function of the central 
baryon density for the sequence as shown in Fig. (5] The max¬ 
imum frequency reached by a non-magnetized and uniformilly 
rotating white dwarf is 1.52 Hz. 


relativistic and magnetic white dwarf is 2.09 Mq, whereas 
the maximum mass achieved by rotation is ^1.45 Mq. 

In Fig. ini we show the relation between the Keple¬ 
rian frequency (/k) and the central density of the star. 
The higher angular velocity, the higher centrifugal forces, 
which pushes the matter outward, and therefore, acting 
against gravity. As a result, the stars are allowed to have 
more mass and, then, increasing the central density. This 
is possible, because the centrifugal forces due to rotation 
(/c oc ril^) have much more effect on the outer layers of 
the star. On the other hand, for non-rotating magnetized 
stars, the Lorenz force acts mainly in the inner layers of 
the star, reducing, and not increasing, the central densi¬ 
ties in these objects as shown in Ref. (l^ . 

Henceforth we will investigate the role played by the 
magnetic field in uniformly rotating white dwarfs . For 
a star with central enthalpy of He = 0.005 c^, whose 
mass is close to the maximum mass in the static case, we 
present results of three different calculations: A) static 
and non-magnetized; B) rotating (with the Keplerian 
frequency) and non-magnetized and C) rotating (with 
the Keplerian frequency) and magnetized. The cases A) 
and B) are presented in the mass-radius diagram in Fig. 
4 andjSJ For the case B), the star rotates with its Ke¬ 
plerian frequency of 0.99 Hz. In addition, to construct 
the models (7), we turn on the magnetic field until that 
the limit of numerical convergence is reached. The re¬ 
sulting magnetic field lines and the electric isopotential 
lines At = const are depicted in Fig. [T] 

In order to study how the Keplerian frequency changes 
with the magnetic field, we perform a calculation for 
different current functions / = /o = const, from zero 
(case B) to the maximum value of the magnetic field as 
shown in Fig. [71 As a result, the Keplerian frequency in¬ 
creases with the magnetic field as shown in Fig. [51 In this 
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FIG. 8. Keplerian frequency as function of the central mag¬ 
netic field for non-magnetized (case B, i.e., /k = 0.99 Hz) 
and magnetized (case C, i.e., Jk = 1.13 Hz) white dwarfs. 
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FIG. 9. Gravitational mass as a function of the central mag¬ 
netic field for a star with central enthalpy He = 0.005 (? for 
the cases B) and C), which are the same stars as shown in 

Fig. El 



FIG. 7. Magnetic field lines (up) and electromagnetic poten¬ 
tial lines (down) for a star with He ~ 0.005 (? and gravia- 
tional mass of -^1.57 M©, with a Keplerian frequency of 1.13 
Hz. The dipole magnetic moment reaches 6.93x 10^® Am^ and 
the magnetic field intensity at the center of this white dwarf 
is 1.87 xl0^'‘G. 


can reach a maximum Keplerian frequency of 1.52 Hz. 
However, in the magnetic case (Fig. |8]), the maximum 
Keplerian frequency is reduced to 1.13 Hz, which corre¬ 
sponds to a white dwarf with gravitational mass of ^ 
1.57 Mq and a central magnetic field of 1.87 xlO^'^G. 


way, equilibrium configurations are obtained for higher 
centrifugal forces and, therefore, if the star can rotate 
faster, it can support higher masses, which explains the 
behaviour observed in Fig. [HI 

According to Figure [HI non-magnetized white dwarfs 


V. CONCLUSIONS 

We computed perfect-fluid magnetized white dwarfs 
in general relativity by solving the coupled Einstein- 
Maxwell equations. We have applied a formalism that 
was developed for neutron stars to rotating magnetized 
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white dwarfs. In our case, the equilibrium solutions are 
axisymmetric and stationary, with white dwarfs endowed 
with a strong poloidal magnetic field. 

The observation of super-luminous la supernovae sug¬ 
gests that their progenitors are super-Chandrasekhar 
white dwarfs, whose masses are higher than 1.4 M©. The 
increasing in the mass may be a result from ultra-strong 
magnetic field inside the white dwarfs. The relevance 
of magnetic fields in enhancing the maximum mass of a 
white dwarf were studied and the results were obtained in 
a fully relativistic framework. We have shown that white 
dwarfs masses increase up to 2.09 Mq for a maximum 
magnetic field strength of ^ 10^® G in the stellar center. 
Thus, magnetic fields can, potentially, be responsible for 
super massive white dwarfs. 

The structure of relativistic, axisymmetric and uni¬ 
formly rotating magnetized white dwarfs were investi¬ 
gated self-consistently and all effects of electromagnetic 
field on the star equilibrium were taken into account. As 
usual, non-magnetized configurations at the mass shed¬ 
ding limit support higher masses as their static counter¬ 
parts. However, we have seen that the magnetic field 
is much more efficient in increasing the mass of WDs 
than rotation. For example, the Keplerian sequence has a 
maximum mass of 1.45 Mq, whereas for the maximum 
magnetic field configuration achieved in this calculation, 
the maximum mass of relativistic WDs is 2.09 Mq, i.e., 
33 % larger than the Chandrasekhar limit, and 2.13 Mq 
in the Newtonian framework. 

We have also shown the increasing of the Keplerian 
frequency (fx) with the magnetic field. The higher the 
magnetic field, the higher the Lorenz force, which, in 
turn, helps the star to support more mass than in the 
non-magnetized case. As a result, if these stars are more 
massive, they can rotate faster (see Fig.[5]and Fig.lH]). In 
this case, the maximum mass obtained for rotating mag¬ 
netized white dwarfs is of ^1.57 Mq, with a Keplerian 
frequence of 1.13 Hz. 


It is to be noted that purely poloidal or purely toroidal 
magnetic field configurations undergo intrinsic instabili¬ 
ties as suggested years ago in Refs, [ssl - ls^ . The na¬ 
ture of this instability was confirmed both in Newtonian 
numerical simulations as in Refs. 5^62| and in general 
relativity framework as in Refs. Analytical and 

numerical calculations have also shown that stable equi¬ 
librium configurations are obtained for magnetic fields 
composed not only by a poloidal component, which ex¬ 
tends throughout the star and to the exterior, but also 
a toroidal one, which is confined inside the star 

In addition, the magnetic field might decay due to 
ohmic effects and, therefore, changing its strength and 
distribution in the star (7l] |. Although we model magne¬ 
tized white dwarfs with a purely poloidal magnetic field 
components, which are not the most general case, we 
have shown, in a fully general relativity way, that these 
stars can considerably increase their masses due to mag¬ 
netic field effects and, therefore, contributing to super- 
luminous SNIa. In a future work, in order to have a more 
complete description of the stars presented here, we in¬ 
tend to take into account magnetic field configurations 
with both poloidal and toroidal components, a B-field 
dependent equation of state and effects of the anomalous 
magnetic moment, too. 
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